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ABSTRACT 

A new model for source counts from 8-1100 /im is presented, which agrees 
well with source-count data and the observed background spectrum. The 
model is similar to that of Rowan-Robinson (2001), but with different evo- 
lution for each of the four assumed infrared template types. The evolution 
is modified in two ways; (i) the exponential factor is modified so that it 
tends to a constant value at late times, (ii) the power-law factor is mod- 
ified so that it tends to zero at redshift Zf, rather than as assumed by 
Rowan-Robinson (2001). I find strong evidence from the 850 and 1100 /im 
counts, and from the infrared background, that Zf — 4-5, with some pref- 
erence for a value at the low end of the range, implying that star-forming 
galaxies at z > 5 are not significant infrared emitters, presumably due to 
a low opacity in dust at these early epochs. 

The model involves zero or even negative evolution for starbursts and 
AGN at low redshifts (<0.2), suggesting that the era of major mergers 
and strong galaxy-galaxy interactions is over. 

Key words: infrared: galaxies - galaxies: evolution - star:formation - 
galaxies: starburst - cosmology: observations 



1 INTRODUCTION 

Source-counts at infrared and submillimetre wave- 
lengths, combined with the spectrum of the infrared 
background, give us important constraints on the 
star-formation history of the universe. First indica- 
tions of strong evolution in the properties of star- 
forming galaxies came from IRAS 60 /im counts 
(Hacking and Houck 1987, Lonsdale et al 1990, 
Rowan-Robinson et al 1991, Franceschini et al 1991, 
1994, Pearson and Rowan-Robinson 1996) 

ISO gave us deep counts at 15 providing 
strong evidence for evolution (Oliver et al 1997, 
Rowan-Robinson et al 1997, Guiderdoni et al 1998, 
Aussel et al 1999, Elbaz et al 1999, 2002, Serjeant et 
al 2000, Gruppioni et al 2002, Lagache et al 2003) and 
useful counts at 90 and 175 fim (Kawara et al 1998, 
Dole et al 2001, Efstathiou et al 2000a, Heraudeau et 
al 2004). 

850 /im counts with SCUBA on JCMT (Smail et 
al 1997, Hughes et al 1998, Eales et al 1999, 2000, 
Barger et al 1999, Blain et al 1999, Fox et al 2002, 
Scott et al 2002, Smail et al 2002, Cowie et al 2002, 



Webb et al 2003, Borys et al 2004, Coppin et al 2006) 
have given important insight into the role of cool dust 
and also strong constraints on the high-redshift evo- 
lution of hyperluminous infrared galaxies. Evidence 
for luminosity evolution of submillimetre galaxies was 
given by Ivison et al (2002). 1200 /im counts with 
MAMBO have been reported by Greve et al (2004) 
and modelled in terms of a strongly evolving popu- 
lation. Recently the AZTEC collaboration have re- 
ported differential counts at 1100 /im (Perera et al 
2008, Austermann et al 2008), which provide even 
stronger constraints on high redshift evolution. 

The detection of the infrared background with 
COBE (Puget et al 1996, Fixsen et al 1998, Hauser 
et al 1998, Lagache et al 1999) provided an important 
constraint on evolutionary models (Guiderdoni et al 

1997, 1998, Franceschini et al 1998, 2001, Dwek et al 

1998, Blain et al 1999, Gispert et al 2000, Dole et al 
2001, Rowan-Robinson 2001, Chary and Elbaz 2001, 
Elbaz et al 2002, King and RR 2003, Balland et al 
2003, Xu et al 2003, Lagache et al 2003). Dole et al 
(2006) used stacking analysis on deep Spitzer counts 
at 24, 90 and 160 fim to estimate the integrated back- 
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ground radiation from sources at these wavelengths. 
With Spitzer we also have deep counts at 8, 24, 70, 
160 /ttm (Fazio et al 2004, Chary et al 2004, Mar- 
leau et al 2004, Papovich et al 2004, Dole et al 2004, 
Le Floch et al 2005, Frayer et al 2006, Shupe et al 
2008), with 24 /im providing an especially complete 
picture of the contribution of individual sources to the 
infrared background. While some pre-Spitzer models 
were quite successful at 70 and 160 fim, none cap- 
tured full details of the 24 /im counts. Lagache et al 
(2004) provided early revised models in the light of 
the Spitzer data. 

Rowan-Robinson (2001) modeled infrared source- 
counts and background in terms of four types of in- 
frared galaxy: quiescent galaxies in which the infrared 
radiation (infrared 'cirrus') is emission from interstel- 
lar dust illuminated by the general stellar radiation 
field, starbursts for which M82 is the prototype, more 
extreme, higher optical depth starbursts with Arp 220 
as prototype, and AGN dust tori. The same evolution 
history, intended to reflect the global star-formation 
history, was used for each galaxy type. This model 
was consistent with counts, luminosity functions and 
colour- luminosity relations from IRAS, counts from 
ISO and available integral 850 pm counts. However 
the advent of deep source-count data from Spitzer 
showed that this model, along with other pre-Spitzer 
models, failed, especially at 24 pirn. 

In this paper I show how the Rowan-Robinson 
(2001) model has to be modified to achieve consis- 
tency with Spitzer and other modern data and give 
predictions for the Herschel and Planck wave-bands. 
Some preliminary results were given in Shupe et al 
(2008). 

A cosmological model with A = 0.7, h —Q.72 has 
been used throughout. 



2 METHODOLOGY 

The philosophy is similar to that of Rowan- Robinson 
(2001), to find a simple analytic form for the evo- 
lutionary function, without discontinuities and in- 
volving the minimum number of parameters. I re- 
tain the same four basic infrared templates: (1) qui- 
escent galaxies, radiating in the infrared through 
emission by interstellar dust of absorbed starlight 
(' infrared cirrus'), (2) starbursts, with prototype 
M82, (3) extreme starbursts with much higher dust 
optical depth, with prototype Arp 220, (4) AGN 
dust tori. These templates, which have been de- 
rived through radiative transfer calculations (Efs- 
tathiou and Rowan-Robinson 1995, 2003, Rowan- 
Robinson 1995, Efstathiou et al 2000b), have proved 
extremely effective in modeling the spectral en- 
ergy distributions (SEDs) of infrared galaxies in the 
ISO-ELAIS (Rowan- Robinson et al 2004), Spitzer- 
SWIRE (Rowan-Robinson et al 2005) and SHADES 
(Clements et al 2008) surveys . The luminosity func- 



Table 1. Evolutionary parameters for each component 



component 


P 


Q 


a 


cirrus 


4.0 


10.0 


0.6 


M82 


2.7 


18.0 


0.994 


A220 


2.7 


13.6 


0.92 


AGN 


2.7 


13.0 


0.94 



tions for each component are derived as in Rowan- 
Robinson (2001), from the 60 fira luminosity function 
via a mixture table which is a function of 60 /im lu- 
minosity (Figure 1). These luminosity functions are 
essentially unchanged from Rowan- Robinson (2001), 
except that the Arp 220 component has been as- 
sumed to be more significant in the luminosity range 
Leo = 10 9 — 10 10 Lq, in order to improve the fit to the 
counts at 850 /im (see below). 

The main new feature of the current work is that 
I now allow each component its own evolution func- 
tion, since a single evolution function for all compo- 
nents is not capable of reproducing the latest counts, 
especially at 24 jj,m. 

In Rowan-Robinson (2001) I assumed a pure lu- 
minosity evolution, with L(z) — <fr(z)L(0), and 

<f>(z) = expQ(t - t)/t (t/t ) P (1) 

where to is the current epoch, the exponential 
factor is essentially the Bruzual and Chariot (1993) 
star-formation history and the power law factor en- 
sures that L(z) tends to zero as t tends to 0. Thus 
to/Q has the meaning of the exponential time-scale 
for star-formation. The peak star-formation rate in 
this model occurs at t/to = P/Q- 

However the 24 /im counts show an extended Eu- 
clidean behaviour at bright fluxes and I have therefore 
allowed for the possibility that the evolution function 
tends asymptotically to a constant value at late times. 
I also allow the epoch where the evolution function 
tends to zero at early times to be tf = t(zf) rather 
than zero, where Zf is to be determined from the fits 
to the counts. This option was found to be important 
at submillimetre wavelengths. Thus I assume 

4>(z) = [ao + (l-a Q )expQ(l-y)}(y p -(y f ) p ) (2) 

where y = {t/to), Vf = t(z f )/t . 

The strategy for determining the parameters P, 
Q, ao for each component was as follows. First we al- 
lowed all parameters to vary while seeking a best-fit 
solution to the 24 /nm counts. This fixed the param- 
eters for the cirrus and M82 components. Keeping 
these fixed we set P for the AGN dust tori to have 
the same value as the M82 starburst component and 
then optimized Q and ao for the dust tori to fit the 
24 fim counts, while requiring that the proportion of 
dust tori dominated sources at S(24) = 1 mjy be no 
less than 20 %, based on analysis of the SWIRE sur- 
vey (Rowan- Robinson et al 2005). These parameter 
choices essentially then determined, without further 
adjustment, the counts at 8, 15, 70 and 160 fim and 
the goodness of fit at these wavelengths is a test of 
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Figure 1. Contributions to 60 fim luminosity function of cirrus (solid), M82 (dotted), A220 (dashed), and AGN dust tori 
(long-dashed) infrared components. 



the broad physical validity of the assumed template 
and luminosity function. 

The predicted counts at 850 (and 1100) /im de- 
pend sensitively on the evolutionary parameters for 
Arp 220 starburst component, but these parameters 
are essentially undetermined at 24 /jm. It became 
clear that the proposal of Efstathiou and Rowan- 
Robinson (2003) that cirrus galaxies could make a 
dominant contribution at S(850) = 8 mjy, was not 
consistent with the count model approach adopted 
here, although they do become dominant at S(850) < 
1 mjy in these count models. The evolutionary pa- 
rameters for the A220 component were adjusted to 
give an equal contribution of A220 and M82 star- 
bursts at S(850) = 8 mjy, since SED modelling of 
SHADES sources (Clements et al 2008) suggests both 
templates are important at these fluxes . This also re- 
quired an increase in the luminosity function for A220 



sources at lower luminosities. We are not able to re- 
duce the contribution of M82 starbursts in the submm 
without destroying the fit to the counts at 8-160 /im. 
Evidence that submillimetre sources do not constitute 
a monolithic, Arp-220-like, population has been given 
by Menendez-Delmestre et al (2007), Ibar ct al (2008 
and Pope et al (2008). 

Adopted values of P, Q, a for each of the four 
components are given in Table 1. Figure 2 shows the 
adopted evolution functions for the four infrared com- 
ponents, with values of z/ of 4, 5 and 10 illustrated. 
The model suggests zero or even negative evolution 
for some components at low redshifts and it would 
be interesting to test this with large area far infrared 
surveys. The physical interpretation would be that 
the era of major mergers and strong interactions is 
over, with the remaining star-formation at the present 
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epoch being relatively quiescent or driven by weaker 
interactions. 

I should emphasize that the assumed evolution is 
applied to all infrared luminosities, and I find no need 
to postulate very strong evolution only of luminous 
galaxies as has been proposed in many source-count 
models in the literature. 

The assumption of a pure luminosity evolution 
is perhaps questionable, since we know that galaxies 
grow by mergers, so that we expect some positive den- 
sity evolution with redshift for more massive galaxies 
and negative evolution for lower masses. However the 
evolution we are trying to describe in eqn (2) refers 
only to star-forming, dusty galaxies. If mergers in- 
volve fragments which are not star-forming or have 
little dust, they will not contribute to density evo- 
lution at these wavelengths. It would be desirable in 
future work to relate evolutionary models to semi- 
analytical models for the merger of dark matter halos 
and the subsequent evolution of star-formation and 
dust, but this is beyond the scope of the present pa- 
per. 



3 FITS TO COUNTS 

Here I discuss the fits to counts at 8-1100 /jm. The 
most challenging wavelength is 24 /mi, with a rela- 
tively Euclidean behaviour at bright fluxes, a very 
sharp steepening of the counts between 3 mjy and 
300 /iJy, and a strong flattening of the source-count 
slope fainter than 200 /tJy. Figure 3 shows my fit to 
the Spitzer 24 /jm differential (Euclidean-normalized) 
counts, indicating the contribution of each of the four 
infrared spectral types. Both cirrus and M82 star- 
bursts are important at bright fluxes, but the strong 
steepening of the counts is attributed solely to M82 
starbursts. AGN dust tori make their strongest con- 
tribution at 1-10 mjy. Arp220 starburst make only a 
small contribution at all fluxes. 

The 15 /im predicted differential counts (Fig 4R) 
fitted to ISO data show similar characteristics to the 
24 /jm counts but the interpretation is quite different 
with cirrus galaxies dominant at all fluxes. This re- 
flects the strong difference between the spectral en- 
ergy distributions (SEDs) of cirrus and M82 tem- 
plates in the mid infrared. The 8 /im predicted dif- 
ferential counts are compared to Spitzer data (Fazio 
et al 2006) in Fig 4L. These appear quite different 
from the 15 and 24 /mi counts, because they are far 
more strongly dominated by quiescent galaxies. 

Figure 5 shows the fit to Spitzer counts at 70 and 
160 /mi differential counts. Here we see dominance by 
cirrus galaxies at bright fluxes and M82 starbursts at 
faint fluxes 

Figure 6L shows the fit to SCUBA differential 
counts at 850 /mi (Coppin et al 2006). Cirrus galaxies 
are predicted to dominate at the very highest fluxes, 
though we have no source-count information at these 



fluxes yet. The Planck mission should remedy this. At 
the fainter fluxes sampled by SCUBA surveys I pre- 
dict that Arp220 starbursts should dominate above 
10 mjy and M82 starbursts below this. This is con- 
sistent with the SED modelling of SHADES sources 
by Clements et al (2008), which shows both templates 
being needed to fit the multiwavelength SEDs. There 
is a strong dependence on the redshift at which in- 
frared starbursts are assumed to switch on, Zf, with 
values in the range 4 to 5 favoured by 850 fim counts. 

Figure 6R shows the fit to the GOODS 1100 
differential counts (Perera et al 2008). These are sim- 
ilar to the 850 /im counts in interpretation, with Zf 
= 4-5 again being favoured. This is at first sight sur- 
prising since there is clear evidence now of galaxies at 
redshifts 6-7. However it is consistent with the 850 fim 
spectroscopic redshift compilation of Chapman et al 
(2003, 2005), which shows no redshifts > 4. That work 
was subject to the constraint that the 850 /im sources 
be detected at radio wavelengths and this could af- 
fect the ability to detect high redshift objects, but 
Chapman et al (2005) argue that this is not a major 
selection effect and that no more than 10 % of the 
submillimetre population could be at z > 4. A few 
individual submillimetre galaxies have been found at 
z > 4 (Wang et al 2007, Younger et al 2007, Schin- 
nerer et al 2008), but there does not appear to be 
a large population of z > 4 submillimetre galaxies. 
The negative K-correction operating at submillime- 
tre wavelengths (Franceschini et al 1991, Blain and 
Longair 1993) means that dust-emitting galaxies of 
redshift 7 should have been clearly detectable in 850 
/xm surveys. Presumably the low metallicity at early 
times means star-forming regions are not opaque to 
dust so star-forming galaxies at z = 5-7 emit purely 
in the optical and ultraviolet. 

The models are also consistent with the MAMBO 
1200 /im integral counts of Greve et al (2004), but in- 
tegral counts do not strongly constrain model param- 
eters because of the non-independence of the errors. 
Herschel and Planck submillimetre counts will pro- 
vide much stronger tests of these models. 

Figure 7 shows a composite plot of dif- 
ferential counts at 8-1100 /im, each wavelength 
horizontally displaced for clarity. Predictions for 
all these wavelengths (and others of interest for 
Akari, Herschel and Planck) are available at 
http://astro.ic.ac.uk/ ~mrr/counts 



4 INFRARED BACKGROUND 

The infrared background spectrum provides a strong 
test of the assumed evolution and SEDs. Figure 8 
shows the fit of my best counts model to the infrared 
background spectrum, showing the effect of assuming 
Zf — 4, 5, 10. The 100-850 /tm background measure- 
ments require Zf — 4-5, with 4 giving the best fit. 
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Figure 2. Evolution rate for cirrus (solid), M82 (dotted), A220 (dashed), and AGN dust tori (long-dashed) infrared 
components. Black curves: adopted maximum redshift of 4; red curves: maximum redshift 5; blue curves: maximum redshift 
10. 



So both the background spectrum and the 1100 (im 
counts are both sensitive to 2/ and require Zf = 4-5. 

The fit to the background spectrum from 15- 
850 /im is excellent. In Fig 8 I have shown the rela- 
tive contribution of the different infrared poulations, 
with cirrus galaxies dominating at 8-24 /im, M82 star- 
bursts dominating to 100 jj,m, and then cirrus galax- 
ies again dominating in the submillimetre, by virtue 
of their dominance of the very faint source counts. 
The Arp220, high optical depth, starbursts, which are 
the dominant contribution to the infrared luminosity 
function above 10 12 Lq, contribute only about 10% 
of the submillimetre background. AGN dust tori con- 
tribute about 0.5 % of the mid infrared background 
and a negligible fraction of the submillimetre back- 
ground. 



5 PREDICTED REDSHIFT 
DISTRIBUTIONS 

Detailed predictions of the models are given at a 
wide range of wavelenghs from 8-1100 /im, includ- 
ing the survey wavelengths of Planck and Herschel 
(astro.ic.ac.uk/ ~mrr/counts/). These include the 
redshift distribution for each infrared type at each 
flux. 

Figure 9 illustrates redshift distributions at 200 
/xjy at 24 /im and at 5 mjy at 850 /im. The 850 /im 
distribution agrees rather well with the Chapman et 
al (2003, 2005) measured distribution. The predicted 
24 fim distribution shows far more z > 1 sources than 
in the distribution derived by Rowan- Robinson et al 
(2008) from their SWIRE photometric redshift cata- 
logue, reflecting the rather shallow optical photome- 
try for most of that catalogue. 
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Figure 3. Euclidean normalised differential counts 24 fim. Solid curve: total counts (black: zj = 4, red(upper): Zf = 5) 
; dotted curve: cirrus; short-dashed curve: M82; long-dashed curve: A220; dash-dotted curve: and AGN dust tori. Black 
curves: maximum redshift 4; red curve: maximum redshift 5, blue curve: maximum redshift 10. Filled circles: SWIRE data 
from Shupe et al (2008), triangles: data from Papovich et al (2004). 



6 DISCUSSION 

The model presented here, based on four basic in- 
frared templates is obviously an oversimplification. 
We know that to understand local galaxies we need a 
range of cirrus components (Rowan-Robinson 1992, 
Efstathiou and Rowan- Robinson 2003). From the 
models of Efstathiou et al (2000) the SED of a star- 
burst varies strongly with the age of the starburst so 
the representation as two simple extremes, M82 and 
Arp220 starbursts, is oversimplified. Finally we know 
that we need a range of dust torus models to fit ultra- 
luminous and hyperluminous galaxies (Efstathiou and 
Rowan-Robinson 1995, Rowan-Robinson 2000, Farrah 
et al 2003). 

However the four templates used here do cap- 
ture the main features of the infrared galaxy popula- 



tion. In fitting the SEDs of ISO and Spitzer sources 
(Rowan- Robinson et al 2004, 2005, 2008), we allow 
individual sources to include cirrus, M82 or Arp220 
starburst, and an AGN dust torus, which gives a 
rich range of predicted SEDs. This works for SCUBA 
sources (Clements et al 2008) and for detailed IRS 
spectroscopy of infrared galaxies (Farrah et al 2008). 

The strongest evolution rate is found for M82 
starbursts. Arp220 starburst peak at a slightly later 
redshift (2 rather than 2.5) and at a slightly lower 
amplitude. AGN dust tori, where the evolution re- 
flects the accretion history, also peak at redshift 2, 
but show a less steep decline to the present than M82 
or Arp220 starbursts, suggesting a gradual change in 
the relative efficiency with which gas is converted to 
stars or accreted into a central black hole. 
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Figure 4. Euclidean normalised differential counts at 8 (LH) (data from Fazio ct al 2004) and 15 (RH) (im (ISO data 
from Elbaz et al 1999 (open triangles), Aussel et al 1999 (filled triangles), Gruppioni et al 2002 (open circles); interpolated 
IRAS data (filled circles) from Verma (private communication). Solid curve: total counts; dotted curve: cirrus; short-dashed 
curve: M82; long-dashed curve: A220; dash-dotted curve: AGN dust tori. 




Figure 5. Euclidean normalised differential counts at 70 (LH) and 160 (RH) fim (filled circles: SWIRE Afonso-Luis et 
al (2008 in preparation), triangles: FLS Frayer et al 2006) . Solid curve: total counts; dotted curve: cirrus; short-dashed 
curve: M82; long-dashed curve: A220; dash-dotted curve: AGN dust tori. Black curves: maximum redshift 4; red curve: 
maximum redshift 5. 



Cirrus galaxies show a much more modest evolu- 
tion rate, peaking at redshift 1, and with an ampli- 
tude only a factor of 3 greater than the present epoch. 
However even in relatively quiescent galaxies like our 
own or M81, the star formation activity was greater 
at redshift 1 than at the present day. This would also 
imply a higher interstellar medium dust optical depth 
in spirals at z ~ 1, an effect that may have been de- 
tected (Rowan-Robinson 2003, Rowan-Robinson et al 
2008). 

Lagache et al (2004) modified their 2003 model 



(Lagache et al 2003) for far infrared counts and the 
background spectrum with changes to the luminosity- 
density evolution and to the 2-30 /itm SED of their 
starburst component. Their earlier model involves 
modest density evolution for normal cirrus galaxies 
from z = to 0.4, then constant luminosity-density 
to z = 5, followed by a slow decline to z = 8. For 
starbursts the luminosity-density rises steeply to z — 
0.5, is relatively flat to z = 4, and then declines to 
redshift 8. Lagache et al (2004) also provide success- 
ful fits to counts from 15-850 fj,m and the background 
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log S(850) (Jy) log S(1100) (Jy) 

Figure 6. Euclidean normalised differential counts at 850 (LH) (SHADES data from Coppin et al 2006) and 1100 (RH) 
/im ( GOODS data from Perera et al 2008). Solid curve: total counts; dotted curve: cirrus; short-dashed curve: M82; long- 
dashed curve: A220; dash-dotted curve: AGN dust tori. Black curves: maximum redshift 4; red curve: maximum redshift 
5; blue curve: maximum redshift 10. 



spectrum, with a philosophy rather different from the 
present paper. The evolution function involves several 
discontinuous changes, a feature I have avoided in the 
present work. I see no necessity to modify our well- 
tested SED models for starburst galaxies. Abundant 
Spitzer-IRS data supports our assumed templates (eg 
Farrah et al 2008). 



7 CONCLUSIONS 

I have presented a new model for source counts from 
8-1100 fim, which agrees well with source-count data 
and the observed background spectrum. The model is 
similar to that of Rowan- Robinson (2001), but with 
different evolution for each of the four assumed in- 
frared template types. The evolution is modified in 
two ways; (i) the exponential factor is modified so 
that it tends to a constant value at late times, (ii) 
the power-law factor is modified so that it tends to 
zero at redshift Zf, rather than as assumed by 
Rowan-Robinson (2001). I find strong evidence from 
the 850 and 1100 /im counts and from the infrared 
background that 2/ = 4-5, with some preference for a 
value at the low end of the range, implying that star- 
forming galaxies at z > 5 are not significant infrared 
emitters, presumably due to a low opacity in dust at 
these early epochs. 

The model involves zero or even negative evolu- 
tion for starbursts and AGN at low redshifts (<0.2), 
suggesting that the era of major mergers and strong 
galaxy-galaxy interactions is over. 

Herschel and Planck submillimetre counts will 
provide much stronger tests of these models. 
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Figure 7. Euclidean normalised differential counts at wavelengths from 8-1100 fim, each displaced by 1 dex for clarity. 
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